BOLTZMANN EQUATIONS FOR MIXTURES OF MAXWELL GASES: 
EXACT SOLUTIONS AND POWER LIKE TAILS 

A.V. BOBYLEV(*) AND I. M. GAMBA(**) 

Abstract. We consider the Boltzmann equations for mixtures of Maxwell gases. It 
is shown that in certain limiting case the equations admit self-similar solutions that 
can be constructed in explicit form. More precisely, the solutions have simple explicit 
integral representations. The most interesting solutions have finite energy and power 
like tails. This shows that power like tails can appear not just for granular particles 
(Maxwell models are far from reality in this case), but also in the system of particles 
interacting in accordance with laws of classical mechanics. In addition, non-existence 
of positive self-similar solutions with finite moments of any order is proven for a wide 
class of Maxwell models. 

To Carlo Cercignani, on his 65th Birthday 



1. Introduction 

In this paper we continue the study of self-similar solutions for various physical sys- 
tems described by the Boltzmann equations with Maxwell collision kernels [1-4]. It was 
understood long ago [5] that the key mathematical tool for such equations is the Fourier 
transform in the velocity space. A detailed analytical theory of the classical spatially 
homogeneous Boltzmann equation for Maxwell molecules was mainly completed in the 
80's (see [6] for a review). There were almost no new essential results in this field during 
the 90's, except for some interesting pure mathematical publications, in particular [7,8]. 

Quite unexpectedly, Maxwell models again became the subject of many publications 
in the beginning of the 2000's. The starting point was an idea to use such models for 
descriptions of inelastic (granular) gases. Inelastic Maxwell models were introduced in 
2000 [1] (see also [9] for the one dimensional case). It was clear that all the analytical 
techniques previously developed for classical Maxwell models can be used in the inelastical 
case almost without changes. Many references to papers published in 2000s by physicists 
can be found in the book [10]. 

One interesting result (absent in the elastic case) was the appearance of self-similar 
solutions with power like tails. It was conjectured in [11] and later proved in [3,4], 
that such solutions represent asymptotic states for a wide class of initial data. On the 
other hand, inelastic Maxwell models are just a rough approximation for the inelastic 
hard sphere model and they give usually wrong answers to the question of large velocity 
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asymptotics (a mathematically rigorous study of such asymptotics for a hard sphere 
model can be found in [12].) 

Some new results in the theory of classical (elastic) Boltzmann equation for Maxwell 
molecules were also recently published in [2,3]: Self-similar solutions (two of which 
were found in explicit form) and the proof that such solutions represent a large time 
asymptotics for initial data with infinite energy, clarification of the old Krook-Wu con- 
jecture [13], etc. It is clear that both elastic and inelastic Maxwell models must be studied 
from a unified point of view. 

An interesting question arises in connection with power-like tails for high velocities: 
Is it possible to observe a similar effect (an appearance of power-like tails from initial 
data with exponential tails) in the system of particles interacting according to laws of 
classical mechanics (i.e. without inelasticity assumption)? This is the main question for 
the present paper. 

We shall see below that the answer to this question is probably affirmative. The key 
idea is to consider a mixture of classical Maxwell gases and to find a corresponding 
limiting case for which such behavior can be in principle observed. 

The paper is organized as follows. First we consider the Boltzmann equation for 
Maxwell mixtures and pass to the Fourier representation (Section 2). Then we study 
a binary mixture and show that corresponding equations formally admit a class of self- 
similar solutions (Section 3). In order to simplify the problem we pass to the limit 
that corresponds to a one component gas in the presence of the thermostat with fixed 
temperature T (Section 4). The general problem can be reduced to the case T = (cold 
thermostat). Then we consider the case of the model cross section (pseudo-Maxwell 
molecules with isotropic scattering) and construct a family of exact self-similar solutions 
(Section 5). These solutions are studied in detail in Sections 6 and 7. The new solutions 
have a lot in common with exact solutions from Ref. [2]. They, however, have finite 
energy and therefore are more interesting for applications. 

We did not try to prove neither existence of such solutions in the more general case 
nor to show that they are large time asymptotic states for a wide class of initial data. 
This is done, as a particular case, in our paper [14] jointly with C. Cercignani. 

Instead, we prove in Section 8 a general statement (applied, in particular, to inelastic 
Maxwell models and elastic models in the cold thermostat) concerning non-existence of 
positive self-similar solutions with finite moments of any order. 

Thus it is sufficient in many cases to prove that such positive solution does exist, then 
it always has just a finite number of even integer moments for all values of parameters 
of the equation. 



2. Maxwell mixtures 

We consider a spatially homogeneous mixture of N > 2 Maxwell gases. Each com- 
ponent of the mixture is characterized by the molecular mass m« and the distribution 
function fa = fi(v,t), i = 1, . . . , N, where »6K 3 and t G IR+ denote velocity and time 
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respectively. The distribution functions are normalized in such a way that 

dvfi(v,t) = pi 



I 

U TO 



gij(cos6) = I — J g{cos6), = ^ , i, j = 1, . . . , N , 



M3 

where pi are the number density of the i th component of the mixture. Note that the 
quantities p^, % — 1, . . . , N are preserved in time. 

The interaction between particles is described by the matrix of Maxwell type differen- 
tial cross- sections 

a i:j (\u\,9) = i^j-s^cosfl) , i,j — l,...,N. 

where \u\ is the relative speed of colliding particles, 9 G [0, n] is the scattering angle. 
In case of "true" Maxwell molecules, i.e. particles interacting with potentials 

Uij = — j- , otij > , 

where r > denotes a distance between interacting particles, the following formulas are 
valid [15] 

1/2 

x rriij J " J rrii + rrij 

We shall assume below the same kind of formulas for gij(cos8) with an arbitrary 
function g(cos9) (pseudo-Maxwell particles). 

The corresponding system of Boltzmann equations reads 

^ = E / dv * duj 9ij (^f) [mm) - /<(«)£(«.)] (2.i) 

where the pair (V,t^) are pre-collisional velocities 

v = {rriiV + rrijV* + mj\u\uj){rrii + rrij)^ 1 , 
K = (rriiV + rrijV* — mi\u\uj)(mi + rrij)' 1 ; i, j — 1, . . . , N 

with respect to the post-collisional velocities (t>,t>*). 
The Fourier transform of Eqs. (|2.1|) - ()2.2|) . for 

m 3 



(2.2) 



Vj {k,t) = I dvfjiv,^- 1 ^; k G 

leads to equations [6] 

dt 



r) N 

^£ (2.3) 



where 

s{y>h t Pj)= j d^dij (~^~) [pi&ij'tfwiKjit) - <Pi(Kt)(Pj{o,t)] (2.4) 



5-' 
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where 



, t rriik + mAk\oj , rn,- 

fc+ = — , kZ = 3 - — (k - \k\oj) 

in; + rrij 3 in, ■ in, 



so, 



k = k± + k {j ; i,j = 1,...,N. 



We note that 



(^) 2 = 4 |^| 2 = |fc| 2 1-4 ^ a 

J \i7ii + rrij / J \ {mi + mjY 

where 

1 / fc-uA 

Then, we consider isotropic solutions of Eqs. ()2.3j) - (j2.4|) 

\k\ 2 



(pi(k,t) = (pi 



2rrij 



and obtain 



with 



' 7 = 1 



r 1 r / 1^|2 \ 

S{{p h {pj) = J dsGijis) (pi ^— (1 - Pij8),tJ 



\k\ 2 rrij \ „ /|fc| 2 



and 



AmiTrij 
(mj + mj) 2 

for i, j = 1, . . . , N. We omit the tildes and denote 

\k\ 2 

x = 

zrrii 

in each of Eqs. (|2.5j) . Then the resulting set of equations becomes 

a w r 1 

—<Pi(x,t) = ^2 dsGij(s) [<Pi(x(l - /3ijs))<Pj(xPijs) - <Pi(x)(pj(0)] 
j= i Jo 



where 



0< Aj =4 7 , 2 < 1 , Pij = Pji, P« = l, i,j = l,...,N. 

(rrii + rrij) 2 



Therefore, the most general system of isotropic Fourier transformed Boltzmann equa- 
tions for Maxwell mixtures reads 

d N 

-q£ = ^lijiViii 1 - Pijs)x)(pj{PijS x) - (pi{x)(pj{Q)) , (2.6) 
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where, for any function A(s), s G [0, 1], 

(A(s))=[ A{s)G{s)ds, G(s) = 4tt#(1 - 2s) , 

Jo 

and (2.7) 

fa~ \>i>l „,,//, . . 
^ = \ ' Pv = ' mi i = 1 '>h3 = h'~,N- 

3. Binary mixture 
We consider a special case N = 2 in Eqs. 1)2.6)1 . ()2.7)1 and denote 

mi = M, m 2 = m, <p x = p x ip{x, t) , yj 2 = P 2 ^(x, t) , 

such that <p(0,i) = ijj(0,i) = 1. Then, we obtain 



r, x 1/2 / \ 1/2 



where 



m J \m 12 J 



AmM mM 
P = i — ; — > m i2 



(m + M) 2 ' iZ (m + M) ' 
(<p,ij>) g = [ dsG{s){cp{{l-Ps)x,t)iP{(3sx,t)-tp{x,t)iP(0,t)} 



(3.i; 



following the notation of ()2.7j) . and thus, 

(<p,tp) = (<p,<p)i, (i^,^) = (ip,ip)i- 

We recall the connection of functions ip(x,t) and ij)(x,t) with corresponding solutions 
/i,2(M>*) °f the Boltzmann equations 1)2.1)1 

p l9 (K,t)= f dv h{\vlt)e- ik -\ p rf(K it )= [ dvf 2 (\v\,t)e- ik -\ 

2M J R3 2m J R3 

The usual definition of kinetic temperatures is given by equalities 

Ti = -± I dv\v\ 2 fi{\v\,t), i = l,2. 

Then, one can easily verify that 

T L (t) = -<p'(0,t), T 2 (t) =-0'(O,t). 

where the primes denote derivatives on x. The equilibrium temperatures T eq of the binary 
mixture reads 

Pi Ti + p 2 T 2 A 
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The relaxation process in the binary mixture described by Eqs. (j3.1j) leads to usual 
Maxwell asymptotics states 

<P ^t^oo exp(-T e(? x) , ip ^t^oo exp(-T eg x) . 

By using Eqs. (J3.1|) . one can easily verify (at the formal level) that 

1/2 



Therefore, 



A = ( ^ ) lis) . is) I r/.sG'(.s).s . 



T\(t) = T eq + - e- Ai , T 2 (t) =T eq + - e~ At 

Pi P2 

A = (Tx(0) - T 2 (0)) , A = A ( Pl + p 2 ) . 

Pi +P 2 



(3.2) 



It is easy to see that Eqs. (j3.1|) formally admit the following class of self-similar 
solutions 

<p(x,t) = <5>{xe- At )e- T ^ x , ^(x,t) = V{xe~ At ) e~ T ^ x . 

It is, however, difficult to investigate such solutions (in particular, to prove that cor- 
responding distribution functions are positive) in the most general case. Therefore we 
consider a simplified problem. 

4. Weakly coupled binary mixture 

If the masses M and m are fixed, then Eqs. (j3.1|) contain five positive parameters 
p i ,aij ,j = 1,2. We shall consider below a special limiting case of Eqs. 1)3.1)) (weakly 
interacting gases) such that 

cti2 — > , p 2 — > oo , p 2 \A*i2 = const . (4.1) 
We assume that the other parameters p x , a u and a 2 2 remain constant and denote 

/2a \ 1 ^ 2 

<p(x,t) = 0{x,t) ,ip(x,t) =ip(x,t) , i=pA— j- J t, 

V J (4.2) 

^(^LY /2 = const. 

Then we formally obtain (tildes are omitted below) 

<Pt = (<P,<p) + 0(<p,i(>) fl , (tj>,il>) = 0, (4.3) 

where it is assumed that the functions (p(x,t),ip(x,t) and their time derivatives remain 
finite in the limit ()4.1j) . The limiting temperatures (see Eqs. ()3.2j) ) are given by the 
identities 

T 2 (t) = T 2 (0) = T eq = const., and 2\(t) = T 2 (0) + (T x (0) - T 2 (0)) e - e/3<s) * 
in the notation ()4.2)) (tildes are omitted). 
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Thus 

(V,-0> = o, = 1, ^'(o,t) = -r a (o) ^^(x,t) = e- T2 (°^, 

and we reduce Eqs. (14. 3 J) to the unique equation for <p(x,t) 

^ = (^^) + ^,e-^°)^. (4.4) 

The general case T 2 (0) > can be reduced to the case T 2 (0) = by substitution 

(p(x, t) = (p{x, t) exp(-T 2 (0) x) . 

Eq. ()4.4jl shows that, in the limiting case (|4.1)l . the second component of the mixture 
plays a role of a thermostat with the fixed temperature T 2 (0), moreover, it is enough to 
consider the case T 2 (0) = (cold thermostat). Then Eq. ([4.4)1 . in explicit form, reads 

Msx) p[(l - s)x] + ^[(1 - /3 s)x\ - V {x) [^(0) + 0^(0)]} , (4 g) 

^(0) = V(0) = 1, 

where the argument t of ip(x,t) is omitted. 

We consider below Eq. ()4.5|) assuming that ip(\k\ 2 , 0) is a characteristic function (Fourier 
transform of a probability measure in M 3 ). Then Eq. ()4.5|) describes a homogeneous cool- 
ing process in the system of particles that interact between themselves and with the cold 
thermostat. Though all interactions are elastic (at the microlevel), such system has much 
in common with the gas of inelastic particles considered in [3,4]. It was proved in these 
papers that the general inelastic Maxwell model has the self-similar asymptotics in a cer- 
tain precise sense. We conjecture the same asymptotic property for solutions of Eq. (|4.5|) 
(see [14] for its proof). In this paper we construct some explicit examples of self-similar 
solutions and show that such solutions have power-like tails for large velocities. 




5. Exact solutions in the Fourier-Laplace representation 

We consider Eq. (|4.5j) with (3 = 1 (equal masses m = M = 1) and G(s) = 1 (isotropic 
scattering). Then Eq. ()4.5|) reads 

ds <p[(l - s)x][<p(sx) + 0] - (l + 0)<p(x), (p(0) = l. (5-1) 

From the physical point of view, this is a model for a mixture of two weakly interacting 
gases, which consist of "particles" with identical masses. One can assume that the sorts 
of "particles" differs, say, by color. 

Our goal is to describe a family of self-similar solutions of Eq. (|5.1|) such that 

<p(x, t) = ipixe-^) , ip(x) = l-ax p , x^O, (5.2) 

The parameters 9 > 0, p > and /iGK will be determined later. From now on x denotes 
the self-similar variable. (The notation -^(xe - ^) should not be confused with one for the 
function ip(x,t) from Sections 3 and 4.) 
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Substituting Eq. (|5.2jl into Eq. (|5.1j) we obtain 

1 f x 

Lixip r (x) - (1 + 9)ip(x) + -ip * Oip + 9) = , ?/>i*V> 2 = / dy 4>i{y)4> 2 (x - y) . 

x Jo 

(5.3) 

This equation can be simplified by the use of the Laplace transform similarly to [2] 

w(z) = C(ip)(z) = / ip(x)e~ zx dx , Re z > z . (5.4) 



provided IVK^)! < Aexp(zox), with some positive A and Zq. 
First, we recall properties of the Laplace transform 

C(xip) = —w'(z) , £(x 2 ip) = w"(z) and C(ip') = zw(z) —tp(0) 

and 

= i(z«)- V(0)) = (zw(z))" • 
Then we obtain the following equation for w(z): 

fi(zw)" + (1 + 9)w' + w + 
Next, we denote by 

= zm;(z) = J dxe~ x ip(^—^ , (5.5) 

so that the above equation is transformed to 

fxz 2 u" + (1 + 6)zu + u{u - 1) = . (5.6) 
The next step is to simplify this equation by standard substitutions. We denote 

z = z q , u(z) = u(z) , (5.7) 
and obtain the equation for u{z) (tildes are omitted) 

fiq 2 z 2 u" + q[l + 9 + /i(q - l)]zu + u(u - 1) = . 

Then, setting 

u(z) = z 2 y(z) + B , £ = const., (5.8) 
we obtain the following equation for y(z): 

fxq 2 z 4 y" + z A y 2 + az 3 y' + (3z 2 y + B(B - 1) = , (5.9) 

where the parameters a and f3 are given by the relations 

a = g(5/ig + 1 + 6 - /i) , 

P = 2B-1+ 4/V + 2q(l + 9- fx) . 

Now, the parameters q and B can be chosen in such a way that a — (3 — 0. 
Thus, we take 

1 + - „ 6^g 2 + 1 



BOLTZMANN EQUATIONS FOR MIXTURE OF MAXWELL GASES 



9 



and obtain that Eq. (|5.9jl is reduced to 



M 2 y" + y 2 + B ^ B z , ~ = 0- (5-10) 

This is the simplest standard form, for which Eq. (|5.fij) can be reduced in the general 
case. Finally, equation (J5.1U)) is of the Painleve type if and only if B = or B = 
1. Otherwise, it has moving logarithmic singularities and does not have any "simple" 
nontrivial analytic solutions (this is just a repetition of arguments on a similar equation 
considered in [2].) 

Therefore we consider the two special cases. In both cases 

6 = n-l-5{tq. (5.11) 

Case 1: B = 0, which implies 

6/V = -1, so that y" = 6y 2 . (5.12) 
Case 2: B = 1, which implies 

6fiq 2 = 1, so that y" = -6y 2 . (5.13) 

The general solution of the non-linear ODE y" = by 2 is expressed in terms of Weier- 
strass elliptic function. Similarly to [2], we can show that just the simplest exact solutions, 
namely 

Case 1: y(z) = (c x + z)~ 2 and Case2: y(z) = -(c 2 + z)~ 2 , 

with constant C\ and c 2 , are solutions to Eqs. (|5.12j) and (|5.13|) respectively, that lead to 
solutions of Eq. (|5.6|) satisfying appropriate boundary conditions at infinity. 

Coming back to the original notation (see the transformations (15 .7J) . ()5.8|) ). we obtain, 
therefore, two different exact solutions of Eq. ()5.6|) 

Case 1: u(z) = (1 + c^y 2 for 6fiq 2 = -1 ; 

(5.14) 

Case 2: u{z) = 1 - (1 + c 2 z~ q y 2 for 6fxq 2 = 1 . 

We have that, in both cases, the coupling constant 9 must satisfy (j5.11|) . 

The constants c\ and c 2 are determined by the boundary conditions as follows. The 
given asymptotics in ()5.2|) for ip(x), x — > 0, leads to the asymptotics for u(z), as defined 
in (j5.5|) . at infinity 

/ x , b 

H z ) - 1 + — , z 00 , 

zp 

where b is a non-zero constant. Recalling that |V'( ::c )| — 1 fc> r positive solutions of the 
Boltzmann equation, we can assume, without loss of generality, that 

u(z) ~ 1 , z — > oo . (5.15) 
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Then we obtain, for the above two cases in (|5.14j) . the following formulas satisfying 
the boundary condition (|5.15j) at infinity: 



Case 1: q = p; u(z) = ( 1 + -z p 



-2 



6p 2 







V 



Case 2: q = -- , u(z) = 1 - (1 + z 



p/2x-2 



2 

2 ' " v ' / ~ v " ' " ' 1 r ~ 3p 
The result can be formulated in the following way. 



(3p-l)(l-2p) 
2p 2 

(3p+l)(2-p) 
3p 2 



(5.16) 



Proposition 5.1. Eq. (|5.1jl has exact self-similar solutions (I5.2jl satisfying the condition 



x> J 



r(p + i) ' 

/or £/ie following values of the parameters 6(p) and fi(p) 



x — > , p > 



1: /x(p) 



6p 2 



2: Mp) = ^f, 



0(p) = 
0(p) 



(3p-l)(l-2p) 
6p 2 

(3p+l)(2-p) 
3p 2 



(5.17) 



T/ie solutions of Eq. (15.1)1 are gryen equalities 



i = 1,2 



(5.18) 



wift ^1,2(2) /rom ()5.16|) . /or cases 1 and # respectively. 



The solutions have a physical meaning if, both, 6 > and ^(^j - ) is the Fourier 
transform of a positive function (measure). 

The first condition leads to inequalities ~ < p < ~ in the case 1, and to < p < 2 in 
the case 2. The second condition will be discussed in the next section. 



6. Distribution functions 
First we evaluate the inverse Laplace transforms (|5.18jl . In the case 1 we obtain 



i + z - 

z\ 2 



The general formula from [17] leads to 

^i(x) = $(2- 1/p x) , $(x) = ZT 1 



1 1 
- < p < - 
3~2 



:i+z 



-p\-2 



^smp-K I ~~ -xs- 1 /* (l + scosp7r) 



P n Jo 



(6.1) 



'1 + s 2 + 2s cosp7r) 2 
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(3p-l)(l-2p) 
6p 2 



(6.2) 



in Eq. (|5.1|) . We note that 9 = for p = |, ~. Eq. (|5.1|) in such cases is the Fourier 
transformed Boltzmann equation for one-component gas. The exact solutions ()6.1|) with 
were already obtained in [2]. Eq. (|6.1|) therefore yields a generalization of 



V 



i i 

2' 3 



these solutions to the case of binary mixture (with equal masses) provided the parameter 
9 is given (for given p e [1/3, 1/2]) in Eq. ()6.2|) . The corresponding solutions of the 
Boltzmann equations are positive solutions with infinite energy. All their properties can 
be studied in the same way as in [2]. 

The case 2 is more interesting since it includes also solutions with finite energy. The 
inverse Laplace transform 



C '{z [1 (l + ^)2 ] } 



can be evaluated in the following way. We denote 

1 



[1 + zp/ 2 ) 2 



then 



FX /*00 PC 

1- / dy*(y) = / dy*(y) , / 

J0 Jx JO 



(6.3) 



The function ty(x) can be expressed through the integral 



1 

2~7ri 



dz e a 



(1 + zp/ 2 ) 2 ' 



< p < 2 



where the contour C lies around the negative half of the real axis (see [17] for details). 
Then we obtain 



where 



A{x) = Im [1 + (re-^f' 2 ] 



m( x ) = - f dre- rx A(r) , 
7T Jo 

2 2W 2 (l + r p / 2 cosf )sinf 



( 1 + rP + 2r p/2 cosf ) 2 
Coming back to the function ip2{x) ()6.3)1 . we obtain 



dre rx A(r) 



(6.4) 



The final result 

i> 2 {x) = 



4 sin ^ 
pn 



ds{l 



s cos 



;i + s 2 + 2scosf) 2 



< p < 2 



is obtained by substitution s = r p / 2 in the integral ()6.4|) . 
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We remind to the reader that the corresponding distribution function f(\v\,t), that 
solves the Boltzmann equation, reads 

f(\v\,t) = e 3 ^F(\v\e^ 2 ) , = , (6-5) 

where 

-i k-v 



F[F] = [ dvF{\v\)e- il r r 2 

Noting that 



e~ T ^ =F[M T (\v\)] , M T (|u|) 



(2ttT) 3 / 2 ' 

we obtain the integral representation oi F(\v\): 

This function is obviously positive for < p < 1. We note that the solution (|6.5|) - ([6.6j) 
corresponds to the value 

g (3p + l)(2-p) 
3p 2 

in Eq. (lo~TJ) . 

7. Solutions with finite energy 

We consider in more detail the most interesting (for applications) case p — 1 in 
Eq. (pTTIj). Then 

A* = |, * = /(M,*) = e'F(Me</ 3 ) , (7.1) 



where 



We denote 



4 p , exp(-|t;| 2 /2s 2 ) , . 



7T 



F(M) = (2**)-W*(y) , (7.3) 



where 



r 2 „„, 1 f 00 , r 2 e~ r 







*(ir7F^ = 7 jf *(TTF' (7 ' 4) 



y 

Asymptotic expansion of <&(y) for large positive y follows from integration of the formal 
series 

( 1 + ^gM)> + 1 >Q' 

Thus we obtain from Eq. ()7.4j) 



Hy),^-l) Jn+1)l : + 2)l , V-oo (7.5) 
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In order to describe a behavior of $(y) for small positive y we transform the first 
integral in Eq. (J7.4j) in the following way 



dr 



(r-1) 



-(r-l)y _ y 



-ry 



Noting that 



we obtain 



where 



ey\E (y) - 2E 1 (y) + E 2 {y)\ , E m (y) = / dr 



E (y) = — , E 2 (y) = e~y - yE 1 (y) 



Hy) = l + --(2 + y)e y E 1 (y) 

y 



EM 



ds - — 

s 



ds - h / ds ' 



lay . 



We note that 



oo e - s rl 
ds V I ds 



dse s In s = —7 



where 7 ~ 0, 577 is the Euler constant. Therefore 

E x {y) = -(7 + In j/) + ds 

Jo 

and 

i r 

e^(7 + lay) - 



$(y) = l + - + (2 + y) 

y 



(7.6) 



Thus, the asymptotic equality from ()7.5|) and the formula ()7.6|) describe the behavior of 
the distribution function F(v) in ()7.3|) . for large and small values of \v\. We obtain 



• 2 x5/2 1 

F(\v\) = 2(-) ± 



1 + 



u — * OO 



f(H) = ^ |4[i + 2H a iBH + o(H a )] 



(7.7) 



M -»• . 



All the exact solutions can be generalized to the case of the thermostat with finite 
temperature T (see Eq. (14. 4j) with ^2(0) = T). Then Eq. (|5.1jl is replaced by the following 
equation 

^ = J ds{<p(sx)<p[(l-s)x]-<p(x)cp(0)} + 9 ds {<p[(l - s)x]e- Tsx - <p(x)} , <p(0) 

This equation can be reduced to Eq. (|5.1|) by substitution 

(f(x,t) = 0(x,t)e~ Tx . 
The corresponding self-similar solutions read 

<p(x,t) = ^(xe-^) e - Tx , 
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where ip( x ) satisfies Eq. (jf>.3j) . 



8. Self-similar solutions and power like tails 
We consider in this section a more general class of equations for the function ip{x,t): 

d(f 
~dt 

<p(0) = 1 



dsG(s){<p(a(s)x)ip[(b(s)x]-<p(x)<p(0)} + 8 / ds H(s) {<p[c(s)x] - ip(x)} 



J 



with non-negative functions G(s), H(s), a(s), b(s) and c(s) with s £ [0, 1]. We also assume 
that G(s), H(s) are integrable on [0, 1], and 

a(s) < 1, 6(s)<l, c(s)<l; 0<s<l. 
The function tp(x,t) is understood as the Fourier transform 

n \k\ 2 ,t)= f dvf(\v\,t)e- tk - v ; keR d , d=l,2,..., ( 8 .2) 

of a time dependent probability density f(\v\,t) £ M. d . 

Eq. (|8.1j) allows to consider from a unified point of view two different kind of Maxwell 
models: 

I Inelastic Maxwell models [1,3], where 

H(s) = 0, a(s) = sz 2 , b(s) = l-sz{2-z), z = — ^,0<e<l. (8.3) 

II Maxwell mixtures described by Eq. ()4.5|) . where 

H(s) = 6G(s), a(s) = s, b(s) = 1 - s, c{s) = l-(5s. (8.4) 

The condition of integrability of G(s) and H(s) can be easily weakened. We do not 
do it here in order to simplify proofs. Our main goal in this section is to prove, roughly 
speaking, that self-similar solutions (distribution functions f(\v\,t) from Eq. (J8.2)0 have 
power-like tails. More precisely, we are going to show that such distribution functions 
can not have finite moments of any order. 

Eq. ()8.1|) admits (formally) a class of self-similar solutions 

<p(x,t) = tp{xe-^) , yu>0, (8.5) 
where <f(x,t) satisfies the equation 



— jjxip' = j dsG(s) [ip(a(s)x)ip((b(s)x) — ip(x)} + \ dsH(s) [ip(c(s)x) — ip(x)] 

V>(o) = i, //>o. 



•6) 
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The corresponding function f(\v\,t) G M d (see Eq. f!8.2jl ) reads 

f(\v\,t) = e d fF(\v\ef), V(l*f) = / dv F(\v\)e~ ik - V . (8.7) 



Our goal is to prove the following general fact: if such a function -F(H) > (generalized 
density of a probability measure in M. d ) does exists, then it can not have finite moments 
of all orders. We assume the opposite and represent the integral in Eq. (|8.1|) as a formal 
series 

°° (—-\) n 

m\ 2 ) = J2-pru a ^ m n\ k \ 2n i 



(2n)\ 

n=0 v ' 



where 



m n = dv F(\v\)\v\ 2n , n = 0,l,... ; 

a„(l) = l, a n {d) = _!_/' du{u'-uf n , ifd>2. 

where IS^ -1 ) is the usual measure of the unit sphere in M. d , uj' G S^ 1 is an arbitrary 
unit vector. 

Hence we obtain 



l*| 8 ) = £(-l) n ^ 



\k 



2)i 



n=0 



Til 

■00 = 1, 4>n = j^-^ra n (d) m n , n = l,2, .... 

The convergence of the Taylor series (|8.8jl is irrelevant for our goals. The only impor- 
tant point is that ip(x) is infinitely differentiate for all < x < oo (see any textbook in 
probability theory, for example [16]) and 

V (n) (0)=^n>0, 71 = 0,1,.... (8.9) 

On the other hand, the equations for ip n can be easily obtained by substitution of the 
series flESJ into Eq. (fFTKj) . Then 

$, = 1, 0i[/i-A(l)] = 0, 

n-l 

il) n \jin- \{n)\ + ^ £(£;, n - k)^ k ^) n _ k = , n = 2, 3, . . . , 



fe=i 



where 



A(n)= / dsG(s)[l -a n (s) -b n (s)} + / ds#(s)[l - c n (s)] , 

* 

g( k j)=( k + l \f dsG(s)a k (s)b\s), k,l = \,2,.... 
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Now we can use conditions f)8.9j) . First we obtain /i = A(l) and recall that // > by 
assumption (|8.5j) . Then we note that Q(k, I) > for all k,l — 1, 2, Therefore 

V>n[-wn + A(n)l > ^->/i>0, n = 2,3,.... 

n 

On the other hand, 

A(n) < / ds[G(s) + #(«)] < oo, 

o 

and, therefore, we get a contradiction. 

Thus, the following statement has been proven. 

Proposition 8.1. Equation flESJ), where fi > 0,G(s),H(s) G L+[0,1],0 < a(s) < 1,0 < 
b{s) < 1, < c(s) < 1 ; s G [0, 1], (ices not /iat»e infinitely differentiate at x = solutions 
satisfying conditions (j8.9|) . 

Corollary 8.2. TTie corresponding probability density F(\v\) from (18. 7|) cannot have finite 
moments of all orders. 

We can now apply the result to the inelastic Maxwell models (j8.3J) and conclude that 
the similar statement proved in our first paper [1] on that subject (see [1], Section 5, 
Theorem 5.1 ) is valid not just for almost all, but for all values of the restitution coefficient 
e G [0, 1]. Consequently we can revise now a statement from [3] related to existence of 
self-similar solutions with finite moments of any order for a countable set of values of e 
from the interval set [0,1] (a possible logarithmic singularity was missing in the sketch of 
proof of Theorem 7.2 in [3]). In fact the solution constructed in [3] has a finite number 
of moments for any < e < 1, without exceptions. 

On the other hand, the above Proposition ^. ll can be applied to Maxwell mixtures (j8.4j) . 
It shows that any physical (i.e. with a positive F(\v |) in Eq. ()8.7^ solution of Eq. ()8.6|) 
corresponds to the distribution function -F(H) with a finite number of even integer 
moments. The exact solution constructed in Section 5 can be considered as a typical 
example. 
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